Spin polarization and effective mass: a numerical study 
in disordered two dimensional systems 
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We numerically study the magnetization of small metallic clusters. The magnetic susceptibility 
is enhanced for lower electronic densities due to the stronger influence of electron-electron interac- 
tions. The magnetic susceptibility enhancement stems mainly from an enhancement of the mass 
for commensurate hllings, while for non-commensurate fillings its a result of an enhancement of the 
Lande g factor. The relevance to recent experimental measurements is discussed. 
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Much attention has recently been given to the study of 
interacting electrons in two-dimensional disordered sys- 
tems, motivated by new experimental observations^. The 
conductance of dilute 2D electron systems show fascinat- 
ing temperature and magnetic field dependencies. One of 
the most intriguing behaviors is the strong decrease in the 
critical magnetic field, B c , needed in order to fully spin 
polarize the system at low densities. These low densities 
are characterized by a high ratio of the intra electron 
interaction energy E c and the Fermi energy Ep. This 
ratio is denoted by r s = E c jEp. In the weakly interact- 
ing regime (r s <C 1) the system behaves as independent 
non-interacting electrons with a magnetic susceptibility 
equal to the Pauli susceptibility \ = 9V?b v (where g is 
the Lande factor, fig is the Bohr magneton and v is the 
density of states at the Fermi energy). As the electron 
density is lowered, measurements show an enhancement 
in the susceptibility^i^iSii. Although there is an on- 
going debate whether these measurements support the 
scenario of spontaneous spin polarization 8 , or just an en- 
hanced magnetic susceptibility*^, it is nevertheless gen- 
erally accepted that a large enhancement of the magnetic 
susceptibility occurs over a large region of densities cor- 
responding to r s > 6. 

An interesting question raised recently^ is whether 
to attribute the enhancement of susceptibility to an in- 
crease of the g factor or to an increase of the effective 
mass (the density of states in a 2D system is propor- 
tional to the effective mass). By a novel experimental 
method Pudalov et. al£ were able to measure separately 
X and to as function of r s . Although there are two dif- 
ferent methods to extract to which do not exactly agree, 
nevertheless, it is clear that m increases much quicker as 
functions of r s than g. Recent measurements by Shashkin 
et. al. 8 strengthen the case for a strong dependence of 
to on r s while g turns out to be constant. This result 
is in contrast to theoretical considerations based on the 
Fermi liquid picture which predict a strong enhancement 
in gaLiLii, but has some similarities to Wigner crystal- 
lization scenariosiSiii. 

Since the development of a theoretical description of 
interacting electrons in disorder systems turns out to be 
rather intricate, a lot of effort has been invested in nu- 



merical studies. For example, the influence of electron- 
electron interactions on the persistent current and con- 
ductance has been extensively studied for small metal- 
lic cluster a 14 ! 1 ^ 1 ^ 1 ^ 18 ! 1 ?: 2 ?! 21 . In this paper we will use 
exact digonalization of small clusters in order to investi- 
gate the influence of r s on x, 9 and m of disordered sys- 
tems. The magnetic susceptibility is strongly enhanced 
as electron-electron interaction increases. Nevertheless, 
the susceptibility shows no clear signs of divergence at 
any finite value of r s . While for commensurate fillings 
the enhancement in x is mainly driven by an enhance- 
ment in the mass to, for non-commensurate systems the 
g factor is strongly enhanced by interactions. 

The disordered interacting cluster is represented by the 
following tight-binding Hamiltonian: 
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+U ^ ( n k,r,°- ~ K )( n i,p,cr> - K)s/\?k,j ~ ?1,pI (!) 



is an electron 
, the number 



where f = (fc, j) denotes a lattice site, a\ -. a 
creation operator (with spin a = —5,+^) 

operator is nkj-a = a>k,j;a a k<r,<" e k,j is the site energy, 
chosen randomly between — W/2 and W/2 with uniform 
probability, V — 1 is a constant hopping matrix element, 
K = n is a positive background charge equal to the elec- 
tronic density n = N/M (where N is the number of elec- 
trons and M the number of sites) and s is the lattice 
constant. The electron-electron interaction is composed 
of the on-site Hubbard interaction Un between electrons 
of opposite spin and the and the long range part U. The 
value of U is related to the electronic density and via 
U = V\J\nnr s and we chose C/h = U 2 K Thus, the 
physical content of the density variation performed in 
experiment is captured by controlling the ratio of the 
Fermi energy to the interaction energy achieved simply 
by changing the interaction strengths U . 

We consider systems composed of N — 4 and N = 6 
electrons residing on 4 x 3, 4 x 5 and 6x6 lattices. We 



set W = 8 for the smaller lattices and W = 5 for the 
larger one, so that the single particle properties corre- 
spond to a diffusive system., For each value of the in- 
teraction strength between < r s < 30 the six lowest 
eigenvalues of the sectors with S z = and S z = 1 are 
calculated^ for a hundred different realizations. 

Using the excitation energies it is possible to estimate 
the magnetic susceptibility x an d the density of states at 
the Fermi energy v. Since we are dealing with a cluster of 
finite size for which the levels are discrete, it is reasonable 
to use a discrete definition of the susceptibility 
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where Ei(S z ) is the lowest eigenvalue in the correspond- 
ing S z sector. This is based on the fact that the first 
change in the magnetic moment of the system will occur 
when the applied magnetic field will be strong enough 
to favor the S z = 1 state as a ground state. That will 
occur once g/is-B = Ei(S z = 1) — Ei(S z — 0), resulting 
in M = fi B . A plot of ( X )/(x*) = (Ei(U,S z = 1) - 
Ex (U, S z = 0)> /(E 1 (U = 0, S z = 1) - Eh. (U = 0,S Z = 0)> 
(where (. . .) denotes an average over different realizations 
of disorder) is presented in Fig. ^ It can be seen that x* 
is strongly enhanced as function of interaction. For small 
values of r s the numerical results follow the well known 
Fermi liquid predictions x/x* — 1 — {V2/Tr)r s 23 . For 
higher values of r s , a linear relation X/x* V r s ~ l/ r o 
(see Fig. 0) between the inverse magnetic susceptibility 
and 1/rg is observed. From this relation one might ex- 
trapolate l/x* — at r s — ro, i.e., full spin polarization 
at finite density. Noting that 1/rf corresponds to the 
electronic density n and that 1 /x* is proportional to the 
magnetic field B c needed to fully polarize the system, 
this relation may be rewritten as B c oc n — uq which is 
exactly an empirical relation proposed in Rcf. 8 based on 
experimental observations. Nevertheless, at large values 
of r s the growth in x* seems to tapper off, resulting in no 
clear evidence of divergence in x* for the diffusive regime, 
although partial spin polarization is seen for part of these 
realizations^. 

As in the experimental measurements, the main prob- 
lem in interpreting the enhancement in x* is that it might 
originate from and enhancement in the g* (i.e., the re- 
sult of some exchange mechanism) or it might result from 
a change in the density of states v* which is unrelated 
to the spin. Nevertheless, it is possible to differentiate 
both effects using the following consideration. When a 
strong enough magnetic field is applied on a system, the 
system will switch from the zero magnetic field ground 
state to the first excited state with higher spin. This 
process is depicted in Eq. and involves both the den- 
sity of states and exchange. In the Landau quasi-particle 
description^ X* Ix — m *9* / m 9, where g/g* = (1 + Fff ) 
and m*/m = 1 + Ff/2, and F s ^ are the usual singlet 
(triplet) Fermi liquid parameters. Thus, both the sin- 
glet and triplet Fermi liquid parameters play a role. On 
the other hand, when one excites the system from the 
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FIG. 1: The inverse magnetic susceptibility x/x* as function 
of r 3 for different lattice sizes and number of electrons. The 
line corresponds to 1 — (\/2/iv)r s which is the Fermi liquid 
correction to the susceptibility. As the interaction increases 
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FIG. 2: The inverse magnetic susceptibility x/x* as function 
of l/r 2 s (which corresponds to the electronic density n) for 
strong interactions. The inset focuses on the region of strong 
interactions. The lines correspond to a linear fit of the form 
1/X* <xl/r 2 3 -l/r 2 . 



lowest singlet state to the first excited singlet state, only 
the density of states (i.e., m, or the singlet Fermi liquid 
parameter) plays a role. The latter is accessible via the 
level spacing between the two lowest excitations of the 
S z = sector with total spin zero. They can be easily 
identified since states in the S z = sector which have 
a total spin larger than zero, have the same energy as 
their counterpart in the S z — 1 sector (with no magnetic 
field the energy of an S — 1 , S z — state is equal to the 
energy of a S = 1,S Z = 1 state) . Thus we identify the 
two lowest states in the S z = which are not degenerate 
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FIG. 3: The single electron density of states v* for different 
lattice sizes and number of electrons, filled symbols corre- 
spond to non-commensurate fillings, while empty symbols to 
commensurate ones. Inset: behavior of v* at strong interac- 
tions. The density of states (i.e., the mass) strongly increases 
for the commensurate fillings, while it shows no signs of di- 
vergence for non-commensurate fillings. 



with states in the S z — 1 sector: Ei(S 
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(E 2 (S = 0)-E 1 (S = 0))' 

A plot of v* jv is shown in Fig. |3 while g*/g = v*x/ v X* 
is presented in Fig. 0] A very different behavior be- 
tween the density of states and the g factor is seen. 
There is a strong enhancement of the density of state 
as function of the interaction strength for commensu- 
rate fillings, while the g factor is strongly enhanced for 
the non-commensurate fillings. Thus, for fillings which 
are prone to Wigner crystallization (commensurate fill- 
ings) the mass is enhanced , while for frustrated systems 
(non-commensurate systems) the g factor is enhanced. 
Although both the 4x4 and 6x6 lattices are at a com- 
mensurate filling for N = 4 it is interesting to note that 
the 6x6 system shows a strong enhancement in the ef- 
fective mass beginning at r s ~ 5, while the 4x4 show 
a strong enhancement of m only much latter (r s ~ 30) 
which is in the vicinity of the expected Wigner crystal- 
lization transition^. We shall comment on this later. 

Although it is rather doubtful to claim that these re- 
sults on small clusters are directly applicable for the ex- 
periment performed on macroscopic samples, it is in- 
teresting nevertheless to note that the behavior de- 
picted here is in line with the recent experimental 
observations^^. For not to low electronic densities all 
groups agree that the critical magnetic field needed to 
completely polarize the electronic spin is proportional to 
(n — no), which is similar to the behavior seen in this 
study. Once lower densities are considered there is a dis- 
agreement between the experimentalist. In Ref£ it is 
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FIG. 4: The Lande g factor g* for different lattice sizes and 
number of electrons. The g factor is strongly enhanced for 
the non-commensurate fillings, while it saturates for the com- 
mensurate ones. 



claimed that the linear relation holds even for low densi- 
ties and is a reliable method to find the critical density 
for which full polarization occurs. On the other hand, in 
Refii doubts are raised on the applicability of the linear 
relation to low densities, and whether a critical density 
exist. Our numerical studies of small cluster seem to be 
more in line with the latter, since we see deviations from 
the linear dependence at low densities and do not see 
a clear indication of full polarization at any finite den- 
sity. Nevertheless, our clusters are to small to draw firm 
conclusions on this point. 

Another interesting point is the different behavior of 
the mass and g factor as function of r a between the com- 
mensurate and non-commensurate fillings. If exchange 
interactions lead the to the enhancement of the suscepti- 
bility, the g factor should be strongly enhanced, while if 
Wigner crystallization is responsible for the susceptibil- 
ity behavior the mass should increase. Indeed, the dilute 
commensurate system (6 x 6, N = 4), which turns into a 
Wigner crystal at large r s , shows a strong enhancement 
in the mass much prior to the expected Wigner crystal- 
lization. The non-commensurate fillings (4 x 3, N = 4 
and 4 x 4, N — 6), which do not show Wigner crystalliza- 
tion at large r s due to frustration, show an enhancement 
in g* but not in v* . The commensurate system (4 x 4, 
N = 4) shows a rather mixed behavior, where both g* 
and v* are enhanced. We believe this is due to the high 
filling which leads to a strong influence of exchange cor- 
relations. Nevertheless, once r s is large enough, g* sat- 
urates, while v continues to grow. In the experimental 
systems, which are believed to exhibit Wigner crystal- 
lization, the large enhancement in the magnetic suscep- 
tibility is attributed^ to the enhancement of the mass. 
The samples used in those experiments are in the dilute 
limit and there is no reason for frustration, i.e., they are 
best described by the dilute commensurate system (6 x 6, 
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N = 4). Thus, the numerical results seem to reproduce 
some of the experimental behavior rather well. 

In conclusion, we have numerically studied the mag- 
netization of small metallic clusters. The magnetic sus- 
ceptibility is enhanced by electron-electron interactions, 
i.e., lower electronic densities. This enhancement stems 



mainly from an enhancement of the mass in the commen- 
surate filling. This has some similarities to the experi- 
mental behavior recently measured. 
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